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ABSTRACT 



Aims. We present observations of twelve rotational transitions of H^O, H^O, and H^O toward the massive star-forming region 

NGC 6334 I, carried out with Herschel/HIFI as part of the guaranteed time key program Chemical HErschel Surveys of Star forming 

regions (CHESS). We analyze these observations to obtain insights into physical processes in this region. 

Methods. We identify three main gas components (hot core, cold foreground, and outflow) in NGC 6334 I and derive the physical 

conditions in these components. 

Results. The hot core, identified by the emission in highly excited lines, shows a high excitation temperature of ~ 200 K, whereas 

water in the foreground component is predominantly in the ortho- and para- ground states. The abundance of water varies between 

4 ■ 10~ 5 (outflow) and 10~ 8 (cold foreground gas). This variation is most likely due to the freeze-out of water molecules onto dust 

grains. The H^O/H^O abundance ratio is 3.2, which is consistent with the l8 0/ l7 ratio determined from CO isotopologues. The 

ortho/para ratio in water appears to be relatively low (1.6 + 1) in the cold, quiescent gas, but close to the equilibrium value of three in 

the warmer outflow material (2.5 ± 0.8). 

Key words. ISM: molecules — stars: formation 



1. Introduction 

Water is one of the most important coolants in star-forming re- 
gions, and is thus a key molecule in the process of star forma- 
tion. Water plays a crucial role in the energy balance in protostel- 
lar objects, and therefore its abundance is an important param- 
!* ■ eter, especially in massive star-formation, because it influences 
, ZZ the collapse and fragmentation of cloud cores. Furthermore H2O 
^^ . (either in the gas phase or as ice on dust grains) is thought to 
5_i ' contain a significant fraction of oxygen. Thus well constrained 
water abundances provide key input for astrochemical models of 
star- forming regions. 

Most rotational lines of H2O, including the ground-state 
transition of ortho- and para-lHbO, are not observable from the 
ground. However, previous space missions have already pro- 
vided a glimpse of the water universe (ISO, van Dishoeck & 
Helmich 1996; SWAS, Melnick & Bergin 120051 Odin, Bjerkeli 
et al. |2009t . The derived water abundance varies widely, between 
~ 10~ 4 in warm (> 100 K) gas and ~ 10~ 8 in cold regions (e.g., 
Boonman et al. 2003). These differences in the H2O abundance 



* Herschel is an ESA space observatory with science instruments 
provided by European-led principal Investigator consortia and with im- 
portant participation from NASA 



are caused by the freeze-out of water molecules onto dust grains. 
Unfortunately, all the observations to date suffer from a short- 
age of accessible water lines (SWAS and Odin observed the or- 
tho ground-state transition only, while ISO only detected excited 
lines) and low spatial resolution (> 1'). 

Van der Tak et al. (2010) observed the para ground-state tran- 
sition (In - Ooo) of Hl 6 in the high-mass protostellar object 
DR 21 (Main) with the Herschel Space Observatory (Pilbratt et 
al. 2010). They found water abundances between 10~ 10 and 10~ 7 , 
about a factor 1000 lower than those derived in previous studies. 
This shows that the H2O abundance in star-forming regions is 
not well understood and that studies including multiple lines of 
both ortho- and para-H20 are highly desirable to put better con- 
straints on the water abundance in the ISM. These capabilities 
are provided by the HIFI instrument (de Graauw et al. 2010). 

In this letter, we present high spectral resolution observa- 
tions of twelve H^O, Hj 7 0, and H^O lines in the high-mass 
star-forming region NGC 6334 I, The relatively nearby (1.7 
kpc; Neckel [T9781 l massive star-forming region NGC 6334 har- 
bors sites of various stages of protostellar evolution (Straw 
& Hyland I1989I I. Single-dish continuum observations at sub- 
millimeter wavelength revealed a total mass of 16700 M Q 
(Matthews et al. 120081 1. of which 200 M is associated with 
NGC 6334 I (Sandell 120001 . NGC 6334 I, studied exten- 
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sively over the last decades (e.g., Beuther et al. 120071 [2008; 
Hunter et al. [2006), is a molecular hot core associated with a 
cometary-shaped ultra compact H n (UCHII) region (de Pree et 
al. 1995), which exhibits a spectrum with very many emission 
lines (Schilke et al. 2006, Thorwirth et al. 2003). Furthermore, a 
bipolar outflow (Leurini et al. 2006 , Beuther et al. 2008]) and sev- 
eral masers have been detected (e.g., Kreamer & Jackson 1995, 
Ellingsen et al. 1996, Walsh et al 2007). Observations using the 
SMA showed that the hot core of NGC 6334 I itself consists of 
four compact condensations located within a * 10" diameter re- 
gion, which are emitting about 50% of the continuum detected 
in single-dish observations (Hunter et al. 2006). 



2. Observations 

We observed twelve water lines with frequencies between 
488 GHz and 1.113 THz, which are listed in Table These 
observations were conducted between March 1 and March 23, 
2010, using Herschel/HIFl in the dual beam switch (DBS) mode 
as part of the guaranteed time key program Chemical HErschel 
Spectral Surveys (CHESS). The coordinates of the observed po- 
sition in NGC 6334 I are a 2 ooo : 17 ;, 20 m 53.32' s and fcooo : 
-35°46'58.5". We used the Wide Band Spectrometer (WBS) 
with a resolution of 1 . 1 MHz across a 4 GHz IF bandwidth. The 
spectra shown here are equally weighted averages of the H and V 
polarization, reduced with the HIPE pipeline (Ott 2010), version 
2.6. We exported the resulting Level 2 double side band (DSB) 
spectra to the FITS format for a subsequent data reduction and 
analysis with the IRAM GILDAS package. 

The spectral scans of NGC 6334 I consist of DSB spectra 
with a redundancy of eight, i.e., the lower and upper side band 
are observed with eight different local oscillator settings. This 
observing mode allows the deconvolution and isolation of the 
single sideband (SSB) spectra (Comito & Schilke 2002), which 
we present here. 

The HIFI beam size at the frequencies observed is given 
in Table [2] For the main beam efficiency we assumed a value 
of 68%. The velocity calibration in HIFI data processed with 
HIPE 2.6 is subject to uncertainties of up to ~ 2 km s due to 
an approximate model for the spacecraft velocity used. A com- 
parison with a more accurate spacecraft velocity model indicates 
that the velocity scale for the data presented here has an uncer- 
tainty of < 0.3 km s . 



3. Results and discussion 

The spectra of the twelve observed water lines are displayed in 
Fig. [J The shapes of these lines are in general very complex, 
indicating the complex structure of the source itself. They can 
be divided into three main components: the hot core, the cold 
foreground component, and the outflow. Although the spectrum 
of NGC 6334 I is in general very rich in emission lines, the wa- 
ter spectra presented here are relatively clean. Emission lines of 
other molecules interfering with the H2O lines are only seen near 
H 18 2„ - 2 02 (H3C-O-CH3) and H 18 lj - lot (HNC and 
H3COH). 

Below we discuss the three individual components present in 
the water spectra. 



-6.3+0.3, -0.3 +0.12, and +6.2 + 0.15 kms -1 . All three compo- 
nents appear in absorption against the strong continuum emitted 
by the warm dust in the hot core of NGC 6334 I. The parame- 
ters of the absorption features seen in the 1 1 1 -Ooo and the 1 1 1- lot 
spectra are listed in Table [3] The velocity of the different com- 
ponents are derived from Gaussian fits. However, Gaussian fits 
do not correctly match the depth of the absorption. Therefore the 
other parameters are determined by visually inspecting the spec- 
tra and integrating over the relevant velocity ranges. To derive 
the continuum level we applied a linear baseline fit in the vicin- 
ity of the absorption line. Thus the background radiation con- 
tains not only the dust continuum, but also water emission from 
all the components in the background (hot core and outflow). 
The velocity-averaged optical depths of the individual compo- 
nents derived from the line/continuum ratio of the 1 1 1 - Ooo line 
(ground-state para transition) are 2.1 + 0.83, 1.27 ± 0.29 and 
1.39 ± 0.63 for the -6.3 kms 1 , -0.3 kms" 1 , and the +6.2 kms" 1 
component, respectively. The corresponding optical depths of 
the ground-state ortho transition are 1 .79 + 0.53, 0.94 + 0.20, and 
1.18 + 0.38. The velocity ranges over which the optical depths 
were averaged are given in Table Q] Deriving the optical depth 
from the line-to-continuum ratio is based on the assumption that 
the excitation temperature is negligible with respect to the con- 
tinuum brightness temperature, i.e., that the foreground compo- 
nent is cold with respect to the background. The large uncer- 
tainties arise because the optical depth of deep absorption lines 
is very sensitive to the actual continuum level, which is deter- 
mined with an accuracy of ~ 5% in our HIFI spectra. In ad- 
dition, the frequency of the H2O ln-loi (557 GHz) is close to 
the edges of the bands la and lb. Thus the uncertainty is en- 
hanced by variations in the sideband ratio of the receiver. For 
the present calculation we used a sideband imbalance of 5%, de- 
rived from the differences of the spectra taken in bands la and 
lb. The lack of absorption of H 18 1 1 1 - Ooo and H 18 l i - 1 10 
at +6.2 kms 1 gives lcr upper limits for the optical depth of 
the para- and ortho-H20 ground-state transitions of 4.5 and 6.2, 
respectively, assuming an Hi 6 0/H^ 8 ratio of 500, which con- 
strains the optical depth at the line center even more. 

Because almost the complete continuum signal is absorbed at 
-6.3 and 6.2 km s , the absorbing material must be cold. Upper 
limits for the excitation temperature are 5.5 K and 9 K for the 
0-H2O and P-H2O, respectively. These upper limits are derived 
assuming large optical depths, thus the remaining signal at the 
center of the water lines corresponds to the source function of 
the transition, and therefore reflects the water excitation temper- 
ature. 

In the following analysis we assume that all water molecules 
are in the ortho- and para-ground states. This assumption is justi- 
fied given the very high critical densities of these two transitions 
(2 • 10 8 crrT 3 and 1.5 ■ 10 7 ctrT 3 for the ground-state para- and 
ortho-transitions, respectively). The column densities per veloc- 
ity interval are thus given by the formula 



N = 



87T-T 



v 3 gl 



? 3 A 



3.1. Cold foreground gas 

The foreground material as seen in the ground-state ortho- and 
para-Hl 6 transitions, consists of three velocity components at 



where N is the column density, t the optical depth, v the line 
frequency, and c is the speed of light. A stands for the Einstein A 
coefficient and gi and g u are the degeneracy of the lower and the 
upper level of the transition. Subsequently, we integrated over 
the velocity ranges given in Table [1] to determine the column 
density for each component. For the -6.3 km s _1 component we 
considered only velocities > -6 km s _1 , because at velocities < 
-6 km s _1 the signal may be contaminated by absorption in the 
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Fig. 1. Spectra of twelve water lines in the high-mass star-forming region NGC 6334 I. Top row: Ground-state para transitions of 
H^O, H^O and H 2 v O. Second row: Ground-state ortho transitions of H^O, H^O and H^O. Third row: Water transitions with 
moderate upper state energies. Bottom row: Water lines with high upper state energies (> 400 K). 



outflow. The derived column densities and the ortho/para (o/p) 
ratio of the three components are listed in Table [JJ Because the 
uncertainties of the column densities are dominated by systemat- 
ical errors (mainly by the uncertainty of the continuum level and 
baseline standing waves), the errors given in Table [JJ are maxi- 
mum uncertainties. The o/p ratio seems to be lower than three, 



Table 1. Column densities of ortho and para water in the fore- 
ground material. 



Comp. 


vel. range 
kms~' 


N(o-H 2 0) 
10 12 [cm" 2 ] 


N(p-H 2 0) 

10 12 [cm" 2 ] 


o/p ratio 


-6.3 
-0.3 
+6.2 


-6 to -2 
-2tol 
5 to 8 


31.4 ±9.4 
14.2 ±3.1 
17.8 ±5.8 


19.3 ±7.8 
8.8 ±2.0 
9.6 ±4.3 


1.6 ±0.8 
1.6 ±0.5 
1.8 ±1.0 



the value expected in statistical equilibrium, in all three fore- 
ground components. However, the uncertainties of the derived 
ratios are quite large for these highly saturated lines, and thus 
the derived o/p ratios have to be treated with caution. 

Water is expected to be formed with an o/p ratio of three. 
Subsequently, collisions with atomic and molecular ions (e.g., 
H + , H!) can lead to proton exchanges, and thus thermalization of 
the o/p ratio. In cool, dense gas, like the foreground components, 
an o/p ratio lower than three is thus likely. From the o/p ratio of 
1.6 + 0.5 derived for the -0.3 km s component, which has the 
smallest uncertainty, we calculate a spin temperature of 19+5 K, 



indicating that the temperature of the foreground component is 
~ 20 K. An o/p ratio lower than three has also been measured 
toward Sgr B2(M) (Lis et al. 120101 see this paper for further 
discussion). 



3.2. Hot core 

Spectra of H^O 4 22 - 3 3 i, H* 6 6 24 - In, and H* 6 5 32 - 4 41 , 
the three transitions with upper state energies above 400 K, show 
only a single component at a velocity of about -8 .2 km s -1 , emit- 
ted by the hot core of NGC 6334 I. The parameters of these 
lines are listed in Table |4] Integrated intensities are derived by 
integrating over a velocity range, whereas the other parameters 
are derived by a Gaussian fit. This velocity is about 1.5 km s -1 
lower than the systematic velocity of NGC 6334 I, but is con- 
sistent with the velocity seen in CH / = | - I (van der Wiel et 
al. 120101) and HCO + J - 12 - 11. Interferometric observations 
of NH3 (3,3) and (6,6) reveal a velocity of -8.1 km s -1 for one 
of the embedded cores (mm 2), which shows the highest opti- 
cal depth in highly excited ammonia lines (Beuther et al. 120071 . 
In addition to the high-energy Hj 6 lines the H^O 2\\ - 2o 2 
line (E up = 136.4 K) peaks at -8.6 km s -1 (Tabled, suggesting 
that this optically thin line is predominantly emitted from the hot 
core as well. 

Two of the observed lines, H 2 53 2 - <\\\ and H 2 4 22 — 3u 
are expected to show maser activity over a wide range of physi- 
cal conditions (Maercker et al. 120081 Neufeld & Melnick fTWTT l. 
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therefore we excluded these lines from the subsequent analysis. 
To calculate the excitation temperature and water abundance in 
the hot core, we assumed that the two remaining lines are op- 
tically thin and that the levels are populated according to the 
Boltzmann distribution. Furthermore, we adopted an o/p ratio of 
three, an H^O/H^O of 500, and, based on interferometric ob- 
servations (Hunter et al. 2006) a source size of ~ 10" to correct 
for the difference in the beam size. With these assumptions, we 
derive a water excitation temperature (T ex ) of 217 ± 30 K, and a 
total H* 6 column density of 7.5 ± 1 .0 ■ 10 15 cnr 2 . 

The C l8 Oi = 10-9 line is observed very close to the o- 
H' 5 3 12 - 3 03 line (Av = -202 MHz) and allows us to inde- 
pendently estimate the H2 column density of the hot core. The 
integrated intensity of this C is O line is 8.78 Kkms -1 , which re- 
sults, assuming a similar excitation temperatures for C I8 and 
H2O and local thermal equilibrium, in a C 18 column density 
of 7.2 + 1 • 10 14 cm" 2 . This gives N(H 2 )=3.8 + 0.5 • 10 21 cnr 2 , 
adopting a CO abundance of 9.5 • 10" 5 and a C 16 0/C 18 ratio 
of 500, leading to a water abundance of 2.0 + 0.3 • 10~ 6 . The as- 
sumption of a similar excitation temperature of H2O and C 18 is 
questionable considering the large difference in their dipole mo- 
ments (C 18 is likely tracing a slightly more extended, cooler 
region). Sandell (120001 1 derives a dust temperature of 100 K for 
NGC 6334 I. Assuming 100 K as lower limit for T ex of C 18 
leads to a C 18 column density of 1.6 • 10 15 cm" 2 , and subse- 
quently to a lower limit of 8 • 10" 7 for the water abundance. 



3.3. Outflow 

Outflow features are seen in the ground-state ortho- and 
para-transitions of all three water isotopologues, except the 
H^O 1 10 - lot line, which is contaminated by H3COH and HNC 
lines. Furthermore, the outflow can be seen in the Hj 6 2n — 2o2 
and Hl 6 3 12 - 3o3 lines (P-Cygni profiles). The spectrum of 



blue lobe of the outflow. This fraction is calculated from 



-MB 



S -T] B 



H\ 6 1 10 - lot is ver y interesting because it shows indication of 
two outflow components, a broad pedestal ranging from -90 to 
+60 km s and a narrower, but brighter component from -40 
to +20 km s (see Fig.|2]l. Because only one outflow is seen in 
interferometric observations (Beuther et al. 2008) and the cen- 
tral velocities of the two components are similar, the observed 
spectra likely reflect a temperature or density variation within 
the outflow. 

That the blue outflow lobe shows up in emission in 
H^O lio — lot, but weakly in absorption in H^O 1 10 - loi may 
be explained by a higher excitation temperature of the optically 
thick Hj, 6 line, caused by photon trapping, which in turn leads 
to level thermalization. 

A prominent feature is the absorption against the hot core 
dust continuum in the blue lobe of the outflow seen in all ground- 
state transitions except Hi 6 1 10 - loi ( see Fig-Ell. Because ab- 
sorption features show up even in Hj 7 0, which is » 1500 times 
less abundant than H^O, the H^O lines are likely completely 
optically thick. Thus the signal in the H^O ln-Ooo line at ve- 
locities from -40 km s" 1 to -20 km s" 1 (6.82+0.12 K) is the 
source function times the beam filling factor (S ■ t]b) plus the 
dust continuum (7.72 K; see Fig. [2]) not affected by absorption. 
At velocities higher than -20 km s -1 possible absorption by the 
foreground component makes the results unreliable, and there- 
fore this velocity range was excluded from our analysis. From 
the emission in the red lobe seen in the H^O 1 1 1 - Ooo transition 
we conclude that the S ■ t\b is 0.38 K for Hi 6 0. Therefore only a 
fraction (/=17%) of the total dust continuum is covered by the 



/=1 

where Tmb is the detected signal strength and T c is the level of 
the background continuum. No emission is seen from the red 
lobe of the 1 1 1 -Ooo line of the rarer isotopologues, and therefore 
only an upper limit of 0. 1 K can be given for the corresponding 
S -r\B- Because this upper limit is of the order of the measurement 
uncertainty of the signal of both lines, we calculated the optical 
depth and column density of the Hj 7 1 1 1 - Ooo and H^O 1 u - 
Ooo lines assuming a source function of zero (all molecules are 
in the ortho- or para-ground state). This agrees with the physical 
parameters found in previous investigations (e.g., Tk > 15 K, 
n(H2) > 10 3 cnr 3 ; Leurini et al. 2006). We calculated the optical 
depth for each velocity channel with the formula 



t = — In 



-i\-f)-T c 
f-T c 



and the corresponding column density as described in Sect. 13.11 
The maximum optical depths (vi sr = -25.7 km s _1 ) are 0.70 
and 0.29 for H* 8 and H' 7 In - 00 , respectively. The cor- 
responding column densities are 4.1 ± 0.6 ■ 10 13 cm -2 and 
1.1 ±0.3 ■ 10 13 cnr 2 forp-H^O andp-H' v O. Hence the resulting 
H^O/H^O ratio is 3.7 ± 0.6, which agrees well with the val- 
ues of the 18 0/ 17 ratio reported in the literature (e.g., Wilson 
& Rood 1994). Assuming the same continuum coverage for the 
Hi 7 O 1 10 - loi line, we calculated the o-H^O column density 
in a similar way, which leads to an optical depth of 0.34, which 
in turn resulted in a column density of 2.7 ± 1 ■ 10 13 cirT 2 . Thus 
the o/p ratio in the outflow is 2.5 ± 0.8. The higher o/p ratio in 
the outflow is consistent with the gas being warmer (> 50 K) 
compared to the foreground clouds. Comparing these values to 
the total column density of the outflow determined by Leurini et 
al. (2006), and assuming an H^O/H^O ratio of 500, we derive a 
water abundance of 4.3 • 10~ 5 , a value typical for warm (> 100 K) 
gas. 



4. Summary 

HerschelfHlFI studies of multiple lines of water isotopologues 
give important insights into the physical processes in dense 
molecular material. Our observations of NGC 6334 I indicate 
that 

- H2O lines show complex line profiles, with multiple emis- 
sion and absorption features, originating from multiple spa- 
tial components. 

- the H2O abundance varies between about 4 • 10~ 5 in the out- 
flow, 2 • 10~ 6 in the hot core, and 10~ 8 in the cool foreground 
gas. These abundances are in the range of previously deter- 
mined ISM values. 

- the H^O/H^O and H^O/H^O ratios are comparable to the 
18 0/ 16 and 17 0/ 16 isotopic ratios determined from obser- 
vations of the CO isotopologues. 

- the water o/p ratio appears to be lower than the statistical 
value of three in the cold, quiescent material, whereas it is 
close to three in the warmer outflow gas. The lower o/p ra- 
tio in the foreground component may be explained by proton 
exchange reactions with H + and H3 , which lower the o/p ra- 
tio to a value corresponding to the gas kinetic temperature of 
about 20 K. 
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Observations of many additional water lines will be available 
in the future, when the complete spectral scan of NGC 6334 I is 
carried out. This will allow a more detailed modeling of water 
abundances in this interesting source. 
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Fig. 2. Absorption feature of the blue lobe of the outflow. The spectra of the ground-state transitions of p-H^O, p-H^O, p-H^O, 
and o-Hl 7 are displayed. The red line shows the dust continuum level. 
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Fig. 3. Spectrum of H^ O 1 10-I01- The red line marks the dust continuum level. The foreground components and the two outflow 
components are labeled. The velocity range of the broad outflow (pedestal) is marked in light blue and the range of the narrow 
outflow component is marked in yellow. 
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Table 2. Summary of the observed lines 



Transition 


spin 


Eiow [K] 


y [GHz] 


beam size [" ] 


HfO 6 24 - 7n 


P 


844.2 


488.491 


44 


H 2 O 1 io - loi 





34.3 


547.676 


41 


H 2 O 1 io - lot 





34.3 


552.021 


41 


IT, O 1 io - lot 





34.3 


556.936 


41 


H 2 fi O 5 32 - 4 4 i 





702.0 


620.701 


33 


Hjj 2n - 2 02 


p 


100.7 


745.320 


30 


H 2 O 2n - 2 2 


p 


100.9 


752.033 


30 


Hi 6 4 22 -3 3 , 


p 


410.6 


916.172 


25 


H 2 O 3i2 — 3o3 





196.9 


1097.365 


20 


H- O 1 1 1 - 00 


p 


0.0 


1101.698 


20 


IT O In - Ooo 


p 


0.0 


1107.167 


20 


IT O 1 1 1 - Ooo 


p 


0.0 


1113.343 


20 



Table 3. Parameters of the narrow absorption features seen in the two ground-state H^O transitions. J TmbAv is the intensity 
integrated over the velocity ranges of the individual components (Table [TJ. The continuum level given is an average value for the 
individual components. 



Tine 


Vlsr 


Tmb 


J T MB dv 


cont. level 




km/s 


K 


Kkms" 1 


K 




-6.2+0.3 


0.21+0.17 


6.2+0.31 


7.7+0.4 


1 l l -Ooo 


-0.34±0.1 


1.87+0.17 


7.8+0.39 


7.7+0.4 




+6.1+0.10 


0.33+0.17 


7.8+0.39 


9.0+0.45 




-6.8+1.0 


0.22+0.024 


2.04+0.10 


1.54 + 0.08 


lio-loi 


-0.25+0.1 


0.81+0.024 


2.69+0.13 


1.53 + 0.08 




+6.3+0.10 


0.19+0.024 


2.36+0.76 


1.51+0.08 



Table 4. Hot core line parameters 



Tine 


J T MB dv 


V/,„- 


Av 


Tmb 




Kkms" 1 


kms~' 


kms~' 


K 


H 2 "0 4 22 -3 3 i 


4.96+0.11 


-8.3+0.11 


6.1+0.27 


0.66+0.04 


H 2 6 6 24 -7„ 


0.28+0.04 


-8.1+0.16 


3.8+0.41 


0.084+0.012 


H} 6 5 32 -4 41 


4.80+0.06 


-8.6+0.60 


7.8+0.14 


0.522+0.013 


H, O 2n - 2o 2 


3.2+0.2 


-8.6+0.3 


6.3+0.5 


0.46+0.05 



